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Abstract 

We investigate the role played by the Polyakov loop in the dynamics of the chiral phase transition 
in the framework of the so-called PNJL model in the SU(2) sector. We present the phase diagram 
where the inclusion of the Polyakov loop moves the critical points to higher temperatures, compared 
with the NJL model results. The critical properties of physical observables, such as the baryon 
number susceptibility and the specific heat, are analyzed in the vicinity of the critical end point, 
with special focus on their critical exponents. The results with the PNJL model are closer to 
lattice results and we also recover the universal behavior of the critical exponents of both the 
baryon susceptibility and the specific heat. 

PACS numbers: 11.30.Rd, ll.55.Fv, 14.40.Aq 
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Confinement and chiral symmetry breaking are two of the most important features of 
quantum chromodynamics (QCD). Chiral models like the Nambu-Jona-Lasinio (NJL) model 
have been successful in explaining the dynamics of spontaneous breaking of chiral symmetry 
and its restoration at high temperatures and densities/chemical potentials. Recently, this 
and other types of models, together with an intense experimental activity, are underway to 
construct the phase diagram of QCD. 

Results with two massless quarks in QCD show that, at high temperature, the phase 
transition associated to restoration of chiral symmetry is second order and belongs to the 
universality class of 0(4) spin models in three dimensions [lj]. With small quark masses, 
the second order phase transition is replaced by a smooth crossover, a picture which is 
consistent with lattice simulations 21. Various results from QCD-inspired models indicate 

n rt 

(see e.g. Refs. [3|, |4|) that at low temperatures the transition may be first order for large 
values of the chemical potential. This suggests that the first order transition line may end 
when the temperature increases, the phase diagram thus exhibiting a critical endpoint (CEP) 

nriQ n n 

[5|, 16|, |3| that can be detected via enhanced critical fluctuations in heavy-ion reactions [a, [9| ■ 

10]. 



At the CEP the transition is second order and belongs to the Ising universality class 
In the chiral limit a tricritical point (TCP) is found in the phase diagram, separating the 
second order transition line from the first-order one. 



Recent developments in lattice QCD 



ll| indicate that the CEP is likely to be localized 



by a new generation of experiments with relativistic nuclei (CBM experiment at FAIR), 
suggesting to explore the range of baryon number chemical potential /i# = 100 — 500 MeV. 

In a previous work [7], in the framework of the Nambu-Jona-Lasinio (NJL) model, we 
studied the phase diagram, focusing our attention on the CEP and the physics near it, 
through the behavior of the baryon number susceptibility and the specific heat. 

In this work we study thermodynamic properties of strongly interacting matter using the 
Polyakov-Nambu-Jona-Lasinio (PNJL) model. This extended model, first implemented in 
Ref. 12j, provides a simple framework which couples the chiral and the confinement order 
parameters. The NJL model describes interactions between constituent quarks, giving the 
correct chiral properties; static gluonic degrees of freedom are then introduced in the NJL 
lagrangian through an effective gluon potential in terms of Polyakov loops 3, Q] with 
the aim of taking into account features of both chiral symmetry breaking and deconfinement. 
The coupling of the quarks to the Polyakov loop leads to the reduction of the weight of the 
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quarks degrees of freedom as the critical temperature is approached from above, which is 
interpreted as a manifestation of confinement and is essential to reproduce lattice results. 
We emphasize that the reduction of the weight of the quark degrees of freedom might also 
have an important role for the critical behavior. 

This effect should be more visible in the temperature domain, which can be explained by 
the attractive interactions between the quarks and the effective gluon field which shifts the 
chiral phase transition temperature to high values, allowing for a stronger first order phase 
transition. 

Hence it is demanding to use this improved NJL model to investigate relevant thermo- 
dynamical quantities such as the CEP and the TCP. 

Our main goal is to locate the critical end point in the PNJL model [jjj] and confront the 
results with the NJL one and universality arguments. Based on the fact that the CEP is a 
genuine thermodynamic singularity, being considered a second order critical point, response 
functions like the specific heat and susceptibilities can provide relevant signatures for phase 
transitions. We notice that susceptibilities in general are related to fluctuations through the 
fluctuation dissipation theorem, allowing to observe signals of phase transitions in heavy-ion 
reactions 

HQ- 

The Lagrangian of the SU {2)®SU (2) quark model with explicit chiral symmetry breaking 
where the quarks couple to a (spatially constant) temporal background gauge field (repre- 
sented in term of Polyakov loops) is given by 



Cpnjl = q(i^D fl -m)q+~g s [(qq) 2 + ( q% 7s rq f } - U $[A]; T) . (1) 

The quark fields q = (u, d) are defined in Dirac and color fields, respectively with two flavors, 
Nf = 2 and three colors, N c — 3, and rh = diag(m°, mjj) is the current quark mass matrix. 
The quarks are coupled to the gauge sector via the co variant derivative = <9 M — iA 1 . 
The strong coupling constant gstrong has been absorbed in the definition of A^: A^(x) = 
9 strong A£(x)<y where A% is the SU C (3) gauge field and A a are the Gell-Mann matrices. 
Besides in the Polyakov gauge and at finite temperature A^ = 5qA° = —i5%A . 
The Polyakov loop $ (the order parameter of Z3 symmetric/broken phase transition in pure 
gauge) is the trace of the Polyakov line defined by: $ = jj-((V expi dr A 4 (x,t) )) . 
The pure gauge sector is described by an effective potential U ($[A|, T) chosen to 
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reproduce at the mean-field level the results obtained in lattice calculations: 



T 4 2 6 v 1 4 



To \ /T \ 2 , /T X 



where 

fe 2 (T) = a + ai If J + a 2 ( ^ \ + a 3 . y 

and a = 6.75, ai = -1.95, a 2 = 2.625, a 3 = -7.44, 6 3 = 0.75, b A = 7.5, T = 270 MeV. 

The parameters of the pure N JL sector are fixed at zero temperature as in {f| : the three- 
momentum cutoff used to regularize all the integrals is A = 590 MeV, m u m d = 6 MeV 
and g s A 2 = 2.435. They yield M vac = 400 MeV, m n = 140.2 MeV, f w = 92.6 MeV and 
(uu) 1/3 = (-241.5MeV) 3 . 

Finally with E 2 = p 2 + M 2 the SU{N f = 2) PNJL grand potential is given by fl Q: 

fi(M,M;T,/i) = U ($,$,T) +2g s N f (q iqi ) 2 -2N c N f [ E p 

J A {2ir) 

_ 2N f T [ ( In [1 + &e-&-M T + me~ 2 ^-^ T + e - 3 C B ^V T l 

A (2tt) 3 I L J 

+ In [1 + 3<l>e-( Ei+ M T + 3$e- 2 ( Ei+ M T + e -^ Ei+ ^ T ] J . (3) 

We notice explicitly that at T = the Polyakov loop and the quark sector decouples. 

The baryon number susceptibility and the specific heat are the response of the baryon 
number density p q (T,fi) and the entropy S(T,fi) to an infinitesimal variation of the quark 
chemical potential \i and temperature, given respectively by: 

The baryon number density is given by p q = ^% J p 2 dp (n q (fx, T) — n q (fx, TV) where 



n q (n,T) and n q (p,T) are the occupation numbers modified by the Polyakov loop [19J. 

The PNJL thermodynamic potential is an effective potential depending on three pa- 
rameters: M, <3> and $. These parameters are not independent (nor the corresponding 
phase transitions) since they should verify the mean-field equations dQ/dM = and 
<9f2/<9$ = dVt/d§ = 0. With the two last equations one can compute $ and $ as func- 
tions of M for any value of T and p. Hence we consider that the thermodynamic potential 
is an effective potential depending only on M. 



A common feature shared by NJL and PNJL models is that the thermodynamic poten- 
tial may have two degenerate minima at which two phases have equal pressure and chemical 
potential and can coexist according to the Gibbs criterium. In fact, this pattern is character- 
istic of a first order phase transition: the two minima correspond, respectively, to the phases 
of broken and restored symmetry. The quark condensate can be identified with the order 
parameter whose values allow to distinguish the two coexisting phases. As the temperature 
increases, the first order transition persists up to the CEP. At the CEP the chiral transition 
becomes of second order. For temperatures above the CEP the thermodynamic potential 
has only one minimum and the transition is washed out: a smooth crossover takes place. 

So, while from a qualitative point of view the results of both models are similar, they 
differ quantitatively in two aspects: in the PNJL, the CEP and TCP are pushed to higher 
values of the temperature and the size of the critical region is larger. Let us now analyze 
those results in more detail. 

In the left panel of Fig. [1] we plot the phase diagram for both, PNJL and NJL models. 
In the PNJL (NJL) model the CEP is localized at T CEP = 169.11 (79.92) MeV and fi CEP = 
321.32 (331.72) MeV [p c q EP = 2.76(2.95)p )- In the PNJL (NJL) model the TCP is located 
at T TCP = 207.66 (112.08) MeV and /i TCP = 270.80 (286.05) MeV. We remark that our 



values for the CEP in the PNJL model are closer to the lattice results of [llj and the main 
change, with regards to NJL values, is in T CEP and T TCP , a result which seems natural 
since the effects of the inclusion of the Polyakov loop are expected to be more relevant in 
the domain of the temperature. 

The size of the critical region around the CEP can be found by calculating the baryon 
number susceptibility, the specific heat and their critical behaviors. The size of this critical 
region is important for future searches of the CEP in heavy ion-collisions [l^]. To estimate 
the critical region around the CEP we calculate the dimensionless ratio Xq/x{ ree i where xl ree 
is obtained taking the chiral limit m° = 0. The right panel of Fig. [T] shows a contour plot 
for three fixed ratios {Xq/x{ ree = 2-0, 3.0, 5.0) in the phase diagram around the CEP where 
we notice an elongation of the region where Xg is enhanced, in the direction parallel to the 
first-order transition line. We also observe that the critical region is heavily stretched in the 
direction of the crossover transition line as shown in Fig. [TJ 

The elongation of the critical region in the (T, /i)-plane, along the critical line, is larger 
in the PNJL model (see Xq/x{ ree = 2.0 in the right panel of Fig. [T]). It means that the 
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FIG. 1: Left panel: the phase diagram in the PNJL and NJL models. Right panel: the size of the 
critical region is plotted for Xq/Xq Tee = 2,3,5. 

divergence of the correlation length at the CEP affects the phase diagram quite far from the 
CEP and that a careful analysis including effects beyond the mean field needs to be done 

As seen in 



191 ] (Fig. 7), one of the main effects of the Polyakov loop is to shorten the 
temperature range where the crossover occurs (at fi = the crossover occurs within a range 
of 150 MeV for the NJL model and within 115 MeV for the PNJL one), thus resulting in 
higher baryonic susceptibilities even far from the CEP. This effect is driven by the fact that 
the one- and two-quark Boltzmann factors are controlled by a factor proportional to <3>: at 
small temperature, $ ~ 0, results in a suppression of these contributions. The thermal bath 
being then only produced via the 3-quark Boltzmann factor, our physical interpretation is 
that the bath is colorless, quarks being produced only in triplet necessarily colorless in the 
average because of $ being the order parameter of Z 3 in this effective theory, $ ~ indicates 
a partial restoration of the color symmetry. When the temperature increases, $ goes quickly 
to 1, resulting in a (partial) restoration of the chiral symmetry which occurs in a shorter 
temperature range. In fact, the most striking difference between NJL and PNJL models is 
a faster variation with temperature, around any characteristic critical temperature, of the 
PNJL results. 

The crossover taking place in a smaller temperature range can be interpreted as a 
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TABLE I: Critical exponents (C.E.): the arrow — > • indicates the path in the /x(T)— 

W 

direction to the CEP (TCP) for p < fi CEP (T < T TCP ). 

crossover transition closer to a second order one than in the NJL model. This "faster" 
crossover may also explain the elongation of the critical region compared to the NJL one 
giving raise to a greater correlation length even far from the CEP. 

With this indication of the important role of the entanglement of the chiral and the 
Polyakov loop dynamics on the critical behavior of the QCD phase diagram, it is mandatory 
to investigate the behavior of Xq and C in the vicinity of the CEP and their critical exponents, 
in the framework of the PNJL model. For comparison purposes with the NJL model and 
the universality/mean- field predictions, the calculated critical exponents at CEP and the 
TCP are presented in Table I, and will be discussed in the sequel. 

The phenomenological relevance of fluctuations in the finite temperature and chemical 
potential around the CEP/TCP of QCD has been recognized by several authors. If the 
critical region of the CEP is small, it is expected that most of the fluctuations associated 
with the CEP will come from the mean- field region around the CEP 

In the left panel of Fig. [2j Xq is plotted as a function of fi for three different temperatures 
around the CEP. For temperatures below T CEP we have a first order phase transition and, 
consequently, Xq nas a discontinuity. For T = T CEP the slope of the baryon number density 
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FIG. 2: Left panel: baryon number susceptibility as a function of [i for different temperatures 
around the CEP in PNJL model: T CEP = 169.11 MeV and T = T CEP ± 10 MeV. Right panel: 
specific heat as a function of T for different values of H around the CEP: yF EP = 321.32 MeV and 
fj, = v° EP ± 10 MeV. 

tends to infinity at fi — fi CEP , which implies a diverging \q- F° r temperatures above T CEP , 
in the crossover region, the discontinuity of Xq disappears at the transition line. 

A similar behavior is found for the specific heat as a function of temperature for three 
different chemical potentials around the CEP, as we can see from the right panel of Fig. [2j 

These behaviors of Xq an d C are qualitatively similar to those obtained in the SU (2) 
NJL model Q. As we have already seen, the baryon number susceptibility, Xq-> an d the 
specific heat, C, diverge at T = T CEP and /i = /j cep , respectively 0, 0]- In order to make 
this statement more precise, we will focus on the values of the critical exponents, in our 
case e and a are the critical exponents of Xq an d C, respectively. These critical exponents 
will be determined by finding two directions, temperature-like and magnetic- field-like, in 
the (T — /i)-plane near the CEP, because, as pointed out in 2l|, the form of the divergence 
depends on the route which is chosen to approach the critical end point. 

To study the critical exponents for the baryon number susceptibility (Eq. |4]) we will 
start with a path parallel to the //-axis in the (T, /i)-plane, from lower \x towards the critical 
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321.32 MeV, at fixed temperature T 



CEP 



169.11 MeV. In Fig. [2 we plot x Q as a 
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FIG. 3: Left panel: baryon number susceptibility as a function of | at the fixed T = T . 

Right panel: specific heat as a function of |T — T CEP \ at the fixed fi = fi CEP . 



function of [i close to the CEP. Using a linear logarithmic fit 

lnx 9 = — em \ ji — jJ CEP \ 



Cl, 



(5) 



where the term c\ is independent of /i, we obtain e = 0.66 ± 0.01, which is consistent with 
the mean-field theory prediction e = 2/3. 

We also study the baryon number susceptibility from higher /i towards the critical jj° EF '. 
The logarithmic fit used now is lnx q = — e'ln|/x — fi CEP \ + c[. Our result shows that 
e' = 0.69 ± 0.02 pa e. This means that the size of the region we observe is approximately 
the same independently of the direction we choose for the path parallel to the /x-axis. These 
critical exponents are presented in Table I, where we can see that the critical exponents 
for the baryon number susceptibility are approximately the same for both, PNJL and NJL 
models, and are consistent with the mean-field theory prediction e = 2/3. 

On the other hand, in the chiral limit (where the CEP becomes a TCP), it is found that 
the critical exponent for \ q has the value j q = 0.51 ±0.01, for both, PNJL and NJL models. 
Again, these results are in agreement with the mean-field value (7, = 1/2). 

Now, paying attention to the specific heat around the CEP, we have used a path parallel 
to the T-axis in the (T, /i)-plane from lower (higher) T towards the critical T CEP = 169.11 
MeV at fixed fi CEP = 321.32 MeV. In Fig. [3] (right panel) we plot C as a function of T 
close to the CEP in a logarithmic scale. We see that for the region T < T CEP we have 



a = 0.63 ± 0.02 1 . Contrarily to what happens in the NJL model (see Table I and Refs. |7|), 
this value of a is closer to the one suggested by universality arguments in jo]. 

We also observe, as in in Ref. that, in PNJL (NJL) model, for the region T < T CEP we 
have a slope of data points that changes for values of \T — T CEP \ around 0.3 MeV. We have 
fitted the data for \T-T CEP \ < 0.3 MeV and \T-T CEP \ > 0.3 MeV separately and obtained, 
respectively, the critical exponent a = 0.63 ± 0.02 (a = 0.59 ± 0.01) and ot\ = 0.53 ± 0.01 
(«i = 0.45 ± 0.01), which have a linear behavior for several orders of magnitude (see Table 
I). As pointed out in this change of the exponent can be interpreted as a crossover of 
different universality classes, with the CEP being affected by the TCP. 

In both models, the influence of the TCP is stronger in the specific heat rather than in 
the baryon number susceptibility: the closest distances between the TCP and the CEP in 
both phase diagrams occur in the T-direction ^{j ]TCP — T CEP ) < (fi CEP — fi TCP )). When the 
CEP is approached from above the trivial exponent a' = 0.69 (for both models) is obtained. 

Let us now analyze the behavior of the specific heat around the TCP. As shown in Table 
I, we find a nontrivial critical exponent a = 0.40 ± 0.02 only for the NJL model while for 
the PNJL model a = 0.50 ± 0.01. 

In this work, we have considered an extension of the NJL model which couples chiral and 
confinement-like order parameters. We have found that our model in general reproduces 
important features of the QCD phase diagram as the location of the CEP/TCP. In addition, 
these results confirm the general idea that, in contrast to the NJL model, the PNJL model 
provides a quantitative description of QCD thermodynamics near critical points. In the 
PNJL model, the crossover taking place in a smaller T range can be interpreted as a crossover 
transition closer to a second order one than in the NJL model. This "faster" crossover 
may explain the elongation of the critical region compared to the NJL one giving rise to a 
greater correlation length even far from the CEP. We have also studied the baryon number 
susceptibility and the specific heat around the CEP which are related with event-by-event 
fluctuations of /i or T in heavy- ion collisions. An important observation is that, in the PNJL 
model the obtained critical exponents are consistent with the mean-field values, both for the 
baryon number susceptibility and the specific heat, while for the NJL this is only true for 
the baryonic susceptibility, since for the specific heat a is different from e. 

1 We use the linear logarithmic fit InC = —a In \T — T CEP \ + C2 where the term C2 is independent of T. 
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As the CEP lies in the region expected to be probed by heavy ion experiments, it would 
be interesting to find an experimental signature of such a point. Near the critical point, 
and in particular in the path we choose to study the critical exponents of the specific heat, 
there is a possibility of the spinodal decomposition in the first order phase transition. So, 
the competition between features of the first and second order phase transition in the mixed 
phase can allow for nontrivial effects, such as the above referred ones and to which there 
is no information from heavy-ion collisions. Our numerical results which also includes the 
chemical potential can be relevant to this purpose. 

In conclusion, the results with the PNJL model are closer to lattice results and we also 
recover the universal behavior of the critical exponents of both the baryon susceptibility 
and the specific heat. The PNJL model here discussed, allowing for finite dynamical quark 
masses, can provide a convenient tool to study the QCD phase diagram; it allows to establish 
a convenient link between the lattice results, and the NJL model itself where gluonic degrees 
of freedom are missing. 
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